Spin relaxation in the presence of electron-electron interactions. 
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The D'yakonov-Perel' spin relaxation induced by the spin-orbit interaction is examined in disor- 
dered two-dimensional electron gas. It is shown that, because of the electron-electron interactions 
different spin relaxation rates can be obtained depending on the techniques used to extract them. 
It is demonstrated that the relaxation rate of a spin population is proportional to the spin-diffusion 
constant D s , while the spin-orbit scattering rate controlling the weak- localization corrections is 
proportional to the diffusion constant D, i.e., the conductivity. The two diffusion constants get 
strongly renormalized by the electron-electron interactions, but in different ways. As a result, the 
corresponding relaxation rates are different, with the difference between the two being especially 
strong near a magnetic instability or near the metal-insulator transition. 
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Understanding the various mechanisms controlling the 
rate of spin relaxation is important for the further de- 
velopment of spintronics. In disordered systems the 
dominant spin relaxation mechanism in the presence of 
spin-orbit (SO) interactions is the D'yakonov-Perel' (DP) 
mechanism Due to the spin splitting, A so , induced 
by the SO interaction the electron-spin precesses through 
the angle A so r between collisions with the impurities, 
where r _1 is the momentum relaxation rate. When the 
spin splitting is weak, A so t <C 1, the precession axis 
is changed randomly at each collision. The net preces- 
sion after N such collisions is typically ^/NA so t rather 
than NA so t. Consequently, the time taken for the elec- 
tron spin to randomize is ~ 1/ A 2 so t. The spin relaxation 
is therefore directly related to the transport properties 
of the electrons. Indeed, spin relaxation investigated in 
metallic semiconductor systems clearly show that the DP 
mechanism is the dominant relaxation mechanism in the 
metallic phase, with the relaxation rate decreasingas the 
metal-insulator transition is approached; see Ref. [2fand 
especially the extended discussion of Fig. 3 in Ref. Q. 

It is known that electron-electron (e-e) interactions 
play an essential role in determining both the trans- 
port and the thermodynamic properties near the metal- 
insulator transition in two-dimensional (2d) electron sys- 
tems In particular, the spin susceptibility has been 
shown to behave critically as the metal-insulator transi- 
tion is approached 0, 0, 1 indicating that the e-e inter- 
actions may also have a profound effect on the rate of spin 
relaxation. In this paper, we show that because of the 
electron-electron interactions the spin relaxation caused 
by the DP mechanism gives considerably different rates 
depending on the technique employed. 

Optical orientation methods are used extensively in 
semiconductors to create a spin-polarized population of 
electrons. The subsequent relaxation of the spins is then 



monitored to study both the temporal la, li| and the spa- 
tial dynamics of the spin subsystem [lOj . In the presence 
of impurities, the spin density propagates diffusively at 
low frequencies and large distances. One can find how 
the e-e interactions modify this propagation by studying 
the dynamic (retarded) spin susceptibility: 



*r(g,w) = i / dte iut ([S x (t),S x (0)}) 



(1) 



In the presence of Fermi-liquid corrections it acquires the 
form: (we use Matsubara frequencies) : 



xTia^n) = x°(i + r 2 ) 



Dq 2 + A 2 so r/2 



(l + T 2 )u; n + Dq 2 + A 2 T/2 
D s q 2 + D s (A so /v F ) 2 



uj n + D s q 2 + D s {A so /v F f 



(2) 



Here, \s = + ^2) is the Pauli s pin susceptibility 

enhanced by the Stoner factor (1 + T2) for free elec- 
trons x° is proportional to the single-particle density of 
states v. In the presence of a finite SO interaction spin is 
not conserved (for the reasons discussed above) and cor- 
respondingly Xg X (0, a;) ^ E3- From Eq. J5J, it follows 
that the rate of relaxation of a spin population due to 
the DP mechanism is: 



t s 1 = D s (A so /v F ) 



(3) 



where D s = D/(l + T 2 ) El Q is the spin-diffusion 
constant, and D — v 2 7 t/2 is the diffusion constant of 
2d electrons. (We assume that v F and A so include the 
standard Fermi- liquid corrections (ljj-l 

An alternate way to study the SO scattering time is 
to study the effect of SO interaction on the quantum 
interference correction to the conductivity 0|- The SO 
scattering rate introduces a cutoff l/r so in the triplet 
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part of the interference processes, leading to the well- 
known anti-localization effect 0] ■ I n the case of the DP 
mechanism, the SO rate is: 



t^ 1 = Alr/2 = D(A so /v F ) 



(4) 



Note that, t sc ^ is proportional to D, and not D s as in 
Eq. 0. 

Since the two rates t" 1 and t^ 1 are directly related to 
the behavior of the different transport coefficients, let us 
clarify the relationship between the diffusion constants in 
a disordered electron liquid 0,0]. The charge-diffusion 
constant D c is related to the conductivity a and the dif- 
fusion constant D by the Einstein relation: 



I 2 dn n 



2vD 



(5) 



The spin-diffusion constant D s is in turn related to D by 
the "Einstein" relation for the spin density as: 



Xs D s = X ° S D . 



(6) 



The compressibility dn/dfi controls the static limit of the 
polarization operator just like \s controls the static limit 
of the spin-spin correlation function. Eqs. (0) and © 
taken together reflect the fact that both the charge and 
the spin are carried by the same particles. Because of the 
rather different renormalization of the respective charge 
and spin-diffusion constants, a strong deviation in the 
rates t~ 1 and t~ 1 may occur. The difference between 
the two rates can be especially strong near a magnetic 
instability or near the metal-insulator transition 7] where 
X s — > oo whereby D s — ► and hence T.T 1 — ► 0, while both 
D and remain finite. 

In disordered conductors in 2d the parameters of the 
electron liquid, in particular the diffusion constants D 
and D Sl acquire logarithmically divergent corrections as a 
function of temperature due to the combined action of the 
e-e interaction and disorder . The question arises how 
the different spin relaxation rates are renormalized as a 
result of these corrections. In this paper we show that 
for the practically important case of the Bychkov-Rashba 
SO interaction the expressions observed in Eqs. and 
(0J at the Fermi-liquid level are preserved in the course 
of the logarithmic renormalizations. 

In 2d systems with structure inversion asymmetry, typ- 
ical to heterostructures, the SO interaction has the form 



of the Bychkov-Rashba term : 



2m 



+ a so a ■ (z x p) . 



(7) 



Naively, the SO interaction can be looked upon as a 
momentum dependent "Zeeman" interaction. While the 
Zeeman splitting induced by a magnetic field is strongly 
renormalized in the presence of disorder and e-e interac- 
tions 0, 0] , we find that the momentum dependence of 
the SO interaction radically changes the situation. The 
momentum dependence allows H so to be rewritten (up 
to a constant) in the gauge form: 



(8) 



where, the SO interaction appears as a spin dependent 
vector potential, f = ^(axz), with effective charge p so — 
2ma so . (In terms of p so , the spin splitting at the Fermi- 
surface A so = p so VF-) We show below that p so , unlike 
the Zeeman term, is not renormalized. 

The renormalization of the parameters of the disor- 
dered electron gas is best described by the matrix non- 
linear sigma model H3,|^|. The SO interaction can also 
be succinctly described within this model. The disorder- 
averaged N-replica partition function of the interacting 
problem reads: 



S[Q] 



\) = I DQ e 

7Tf 



d 2 r 



4r 



S[Q] } (9) 
Tr Q 2 - TrlnG" 1 + QTQ , (10) 



where the Q-&e\d is an auxiliary matrix field describing 
the impurity scattering. The Green's function in Eq. I|1(J|) 
is easily deduced from Eq. JSJ: 



1 i 
G- 1 = ie + (i-—{-iV + Pso ff + —Q 
2m 2t 



(11) 



with the last term appearing because of the impurity 
scattering. On expanding the action S[Q] about its sad- 
dle point solution the problem of electron diffusion in 
the field of impurities, e-e interactions, and SO interac- 
tion reduces to the non-linear sigma model: 



S[Q] 



T 



Tr D ( VQ + i Pso [f, Q] - 4zTr (eQ) + QTQ 



(12) 



The matrix Q satisfies the constraints: Q 2 = 1,Q = 



Q\ and Tr Q = 0. The components of Q are defined 
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FIG. 1: Kernels Ku and K v indicating two possible ways to 
insert the spin-orbit vertices (a) into the propagator T>. 




as Qlif"n 2 , where ni,Ti2 are the fermionic energy indices 
with e n — (2n + 1)tvT; i,j are the replica indices and 
a, f3 are the spin indices. The trace is taken over all 
these variables. The energy matrix e = e n 5 nm dijSna. 
The factor z is the frequency renormalization factor |2l| 
that determines the relative scaling of the frequency with 
respect to the length scale; z = 1 for free electrons. 

The fluctuations of the Q-field in the particle-hole 
channel (diffusons) determine the propagators T>(q,u>), 
which in the absence of the SO interaction (p so = 0) 
have a diffusion-like sing ularity T>(q, u) = l/(Dq 2 + zuj). 
These propagators describe the diffusive evolution of the 
charge and spin density fluctuations at large scales. The 
fluctuations responsible for the interference corrections 
are in the particle-particle channel (cooperons). Since 
the Bychkov-Rashba SO interaction is invariant under 
time-reversal, the particle-hole diffuson and the particle- 
particle cooperon propagators are related by charge con- 
jugation |23,|2J|. It therefore suffices to study the renor- 
malization of the spin-orbit scattering in the diffuson 
channel only. 

The spin fluctuations can be classified in terms of the 
total spin, S, of the particle-hole pairs. The interaction 
matrix T = T s and T t represent the multiple scattering 
induced by the e-e interactions in the singlet S = and 
the triplet 5=1 diffuson channels, respectively. They 
are related to the static Fermi-liquid parameters Ti and 
T 2 as T„=Ti- r 2 /2 and T t = -T 2 /2. The presence of 
the quadratic term Dp 2 so (QtQt) in Eq. ((12")) introduces 
in the triplet 5 = 1 diffuson and cooperon propagators a 
cutoff, i.e., a gap, proportional to the SO scattering rate 
t~„ = Dp 2 Q . Note that the linear term proportional to 
Dp so (\7QtQ) in Eq. (fT^|) may be gau ged away by the 
transformation Q(r) — > e lPsoT ' r Q(r) [25|. This will, how- 
ever, generate higher order terms in p so since Pauli ma- 
trices do not commute with each other. The question of 
the renormalization of the quadratic term therefore still 
remains. 

The quadratic term leads to inserting two r-vertices 
into the particle-hole propagator T>. These vertices are 
marked by the solid triangles (A) in Fig. QJ As shown in 
the figure, there are two possible ways indicated by the 
kernels Kh and K v in which this can be done. To deter- 
mine the scaling behavior of the cutoff r s ~ x in the diffuson 
channel we study how both these kernels are renormal- 



FIG. 2: Ti corrections to the K v kernel in Fig. involving 
only one momentum integration. 
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FIG. 3: T2 corrections to the K v kernel in Fig. involving 
only one momentum integration. 



Since each momentum integration involving the diffu- 
sion propagators leads to one power of p = 1 / (2t:) 2 uD in 
2d, the corrections to first order in the disorder strength p 
are evaluated by taking into account diagrams with only 
one momentum. In the following, we demonstrate how 
the scaling equations for r^ 1 are obtained by considering 
the renormalization of the kernel K v as an example. 

Corrections to K v limited to only one momentum inte- 
gration involving the interactions T% and T 2 are shown in 
Figs. and [21 respectively. (Similar corrections exist for 
the Kh kernel as well.) With the exception of the first 
diagram in Fig. |2| these corrections originate from the 
linear SO term iDp so (\j QtQ) applied twice. They lead 
to the renormalization of the cutoff l/r so . The renor- 
malized cutoff 1/t^ can be expressed in terms of the 
renormalized parameters D — > D' and p so — > p' so : 




where SD and Sp so are the corrections to D and p so , 
respectively. Given SD, one can extract 6p so from the 
calculated 1/t' so . 5D corresponds to the well known 
Altshuler-Aronov corrections to the conductivity (lif: 

Then for 5p so , we get: 



4 



Spso 4 


f d 2 q 


f duj 


Pso V J 


(2tt)2 J 


2^ 



V 2 (q,Uj)T 2 - 2V 3 {q,Lu)Dq 2 T 2 + V 4 (q,Uj)ujDq 2 Tl 



(15) 



r 



Note that Sp so depends only on the triplet amplitude. 

The corrections given in Eqs. I|14(l and 1|15|) have so far 
employed the static amplitudes Ti and T 2 . To include 
the effects of dynamical screening, these amplitudes are 
extended via the ladder summation |2l|: U Si t(q,uj) — 
T S:t (Dq 2 + zu>)/(Dq 2 + (z - 2T Stt )ui). The dynamical re- 
summation allows the evaluation of the corrections to 
infinite order in the interaction amplitudes so that p re- 
mains as the only expansion parameter. Remarkably, 
Sp so vanishes after this replacement of the T 2 amplitude 
in Eq. (JTHJl by U t (q,u>): 



Sp BO = 



(16) 



Similar analysis when extended to the kernel Kh and to 
the linear term Dp so iyQfQ) in Eq. (|T^|) leads to the re- 
sult that the gauge form D(\/Q+ip so [f, Q]) 2 of the action 
in Eq. I|12fl with the bare value of p so is preserved under 
renormalization. This implies that the cutoff = Dp 2 Q 
retains the same form as in Eq. |0J with D as the renor- 
malized diffusion constant and with p so unrenormalizcd. 

Note that the diagrams presented in Figs. Inland |3| re- 
semble those that are used for the calculation of the spin 
susceptibility [lii. Il4| with the triangle at the top as the 
starting point and ending at the bottom triangle. This 
analogy is, however, misleading since the SO interaction 
Dp so (S7QtQ) in addition to the Pauli matrices is propor- 
tional to momentum. It is this dependence on momentum 
that makes the calculation of the renormalization of p so 
different from the renormalization of the Zeeman term. 

We now study the renormalization effects on the DP 
spin relaxation rate t s _1 . To this end, we analyze the 
form of the renormalized dynamic spin susceptibility that 
follows from the action defined in Eq. H12|) : 



Dq 2 + Dp 2 so 



(z + T 2 ) uon + Dq 2 + Dp 2 so 



Xs 



D s q 2 + D s p 



0J n + D s q 2 + D s p 2 so 
where, the renormalized spin-diffusion constant: 
D s =D/(z + T 2 ) . 



(17) 



(18) 



Note that Eq. 1|17[) has the same structure as the expres- 
sion for Xs X (ly u ) m the disordered Fermi-liquid with the 
parameter (z + T 2 ) substituted for (1 + T 2 ) in \s and 
D s . It follows from Eq. JUJ that the form of the DP 



To conclude, we have shown that the form of the func- 
tional S[Q] in Eq. (|12fl is preserved under renormalization 
by the combined action of the e-e interaction and disor- 
der. It follows from this observation that the relaxation 
rate of the spin density and that of the cutoff in the S — 1 
cooperon/diffuson propagators are given by Eqs. y| and 
Q , respectively, where the parameters D s and D are the 
renormalized diffusion constants 26] , and the spin split- 
ting p so = A so /vp is unchanged. As a result, under the 
circumstances when Xs diverges, either near a magnetic 
instability (T 2 — > oo), or near the metal-insulator transi- 
tion (z —> oo) 0, the spin-diffusion constant D s and r s _1 
vanish, while the scattering rate r," 1 remains finite. 
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spin relaxation rate, r s 



D s p 2 a , given in Eq. (j3Jl it 



also preserved under renormalization. 
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